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The Auger Engineering Radio Array (AERA) is currently detecting cosmic rays of energies at and above 10 1T eV at 
the Pierre Auger Observatory, by triggering on the radio emission produced in the associated air showers. The radio- 
detection technique must cope with a significant background of man-made radio- frequency interference, but can provide 
information on shower development with a high duty cycle. We discuss our techniques to handle the challenges of self- 
triggered radio detection in a low-power autonomous array, including triggering and filtering algorithms, data acquisition 
design, and communication systems. 
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1. Introduction 

Cosmic-ray air showers emit coherent broadband radio 
pulses in the MHz frequency range, a phenomenon discov- 
ered by Jelley and collaborators in 1965 [1]. Advances in 
O digital processing technology, along with the potential of 
^3 the radio technique to determine cosmic-ray primary com- 
^ position with a high duty cycle [2], have led to renewed 
i_i interest in the technique [3[ . Radio enhancements are now 
planned at several neutrino and cosmic-ray observatories. 

The goals of the Auger Engineering Radio Array (AERA), 
[a 20-km 2 extension of the Pierre Auger Observatory in Ar- 
gentina, are to calibrate the radio signals using hybrid air- 
shower measurements, and to demonstrate the potential of 
( I the technique for future arrays. The first stage of AERA is 
^operational and is successfully measuring air showers of en- 
ergies at and above 10 1T eV in coincidence with the Auger 
^1 surface detector (SD) and fluorescence detector (FD) [j. 
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2. Detector and Data Acquisition Design 



The first stage of AERA, deployed in September 2010 
^ and July 2011, consists of 24 radio detector stations (RDSs) 
deployed on a triangular grid with a spacing of 150 m and 
covering an area of ~ 0.4 km 2 . Each RDS measures the 
radio signal with a dual-polarization log-periodic dipole 
antenna (LPDA), aligned north-south and east- west, and 
sensitive between 27 and 84 MHz. The received signal 
is further amplified and bandpass filtered, and each po- 
larization is digitized at 200 MHz with 12-bit analog-to- 
digital converters (ADCs). High-gain (+49 dB after the 
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Figure 1: The AERA data flow. The GPS timestamps of station-level 
triggers (L2) are sent to a central DAQ, which forms multi-station 
coincidences (L3; shown at left). An L3 request is sent back to the 
stations (not shown), and digitized waveforms of the event are sent 
to the central DAQ for storage (shown at right). 



antenna) and low-gain (+29 dB) versions of each channel 
are recorded to extend the dynamic range. 

A field-programmable gate array (FPGA) in the sta- 
tion electronics forms a trigger when a bandwidth-limited 
pulse is observed in the voltage traces (a "level 1" trig- 
ger). These triggers are passed to a CPU within the RDS, 
which can then optionally apply other trigger conditions 
(the "level 2" trigger). The triggers are timestamped using 
GPS receivers at each RDS. 

The CPU in each RDS communicates via Ethernet over 
a fiber-optic network to a central data acquisition system 
(DAQ). The RDS reports the trigger timestamps to the 
central DAQ. The DAQ software forms the "level 3" trigger 
by using time and spatial coincidences of multiple RDSs, 
and the time-domain waveform data of these events are 
requested from each RDS and stored to disk. Figure [T] 
shows the data flow for each trigger level. 
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3. Triggering and Filtering 

Triggering directly on the radio signal of the air show- 
ers (instead of using particle detectors as a trigger) poses 
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Figure 2: Top: unfiltered AERA polar skymap (2.5h, 40k events). 
The dark circle indicates the horizon. Bottom: AERA cosmic-ray 
events coincident with the Auger SD (48d, 18 events). The direc- 
tions have been smeared with a 3° Gaussian, and the color scale is 
log 10 (event density, a. u.). 

some challenges for the data acquisition, due to man-made 
radio- frequency interference (RFI). The continuous back- 
ground level is set by the radio emission from the Galactic 
plane, but any man-made narrowband transmitters add to 
the level above which one must detect air-shower pulses. 
Additionally, man-made pulsed RFI (from sparking elec- 
trical equipment, airplanes, etc.) can mimic the signals 
from cosmic rays. 

A comparison is shown in Fig. [2] between raw level- 
3-triggered data and self-triggered cosmic-ray events de- 
tected in coincidence with the SD. Since the bandwidth 
and computational resources at each triggering level are 
limited, one of the technical focuses for the first stage of 
the array has been to develop various methods to reject 
RFI in order to minimize efficiency losses from bandwidth 
saturation. 

3.1. Narrowband Filtering 

Before triggering on a radio pulse, it is advantageous 
to increase the signal-to-noise ratio by filtering out any 
narrowband transmitters from the digitized antenna sig- 
nals. This is accomplished in a computationally efficient 
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Figure 3: A recorded cosmic-ray radio pulse (at 1500 ns), before 
IIR narrowband filtering (black) and after (red, pulse offset in time). 
The peak signal-to-noise ratio increases from 2.4 to 7.6. 

manner by using a series of infinite-impulse response (IIR) 
notch filters in the FPGA. The IIR filters operate on the 
time-domain signal, and the output of the filter yi is a lin- 
ear combination of input samples Xj and delayed feedback 
output samples yj from the filter: 

Hi = %i — (2 COSCo>AT • + #i-2 

+ (2rcoscjAr • Vi-i) ~ (r 2 • ^_ 2 ). (1) 

The normalized filter frequency ujn is given by the notch 
frequency fjy and the sampling frequency f$, 

lon = 2tt f N /f s , (2) 

and the width parameter r is a value strictly between and 
1, with higher values giving a narrower response function, 
r = 0.99 is typical for a narrow transmitter. 

A complication in the implementation of the IIR fil- 
ters in high-frequency FPGAs arises from the fact that 
one cannot arbitrarily pipeline the feedback computation. 
We have resolved this by using the scattered look-ahead 
pipelining technique (EJ , which increases the filter complex- 
ity but allows more time for computation in the FPGA. 

The coefficients of the Xj and yj in Eq. [1] can be pre- 
computed for the desired notch frequencies, converted to 
a fixed-point representation, and loaded into the FPGA at 
run-time. The current design allows for four tunable notch 
filters for each polarization. See Fig. [3] for an example of 
the filtering. 

An additional filtering technique we are exploring is the 
use of a frequency-domain median filter. The time-domain 
signal is first converted to the frequency domain with a 
fast-Fourier transform (FFT), and then the amplitude of 
each frequency component fa is replaced by the median 
value in a window of size 2N: fi = median (/^-n, /i+jv)- 
An inverse Fourier transform then converts the signal back 
to the time domain for triggering. The median filter has 
the advantage of removing any number of narrowband 
transmitters, but it is computationally intensive and can 
distort pulse shapes in the time domain. The pre-trigger 
FFT also allows for a number of other strategies for signal- 
to- noise improvement, such as deconvolution of the band- 
pass filter dispersion, or use of a matched filter [6[. 
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Figure 4: The AERA LI trigger scheme. See the text for details. 

3.2. Level 1 Trigger 

Designed to capture transient, isolated radio pulses, 
the digitizer uses algorithms in the FPGA to trigger in 
the time domain on potential signals. At its most basic, 
the trigger is simply a voltage threshold above a baseline. 
However, a number of other parameters are added to this 
trigger to reject RFI. The following trigger conditions are 
applied to select clean, bandwidth-limited pulses: 

1. the voltage rising edge crosses the primary threshold 

2. before the T\ crossing, no other T\ crossings occur 
during the previous time period T prev ; 

3. after the T\ crossing, the signal rising edge crosses a 
secondary threshold T2 , where normally T2 < T\ ; 

4. the number of rising-edge T2 crossings NC within a 
time period T per falls within a specified range NC m i n < 
NC < NC max ; 

5. the time TC between successive T2 crossings is less 
than a maximum value TC max ; 

6. the quotient Q of the pulse maximum P max divided 
by the number of T2 crossings NC falls within the 
range Q min < Q < Q max . 

The thresholds are compared relative to a baseline volt- 
age, determined dynamically with a rolling average over a 
specified number of samples. Samples outside of a specified 
voltage range (i.e. part of a large pulse) are not included 
in the baseline calculation. 

Figure E] shows the trigger parameters graphically. Each 
condition is designed to reject man-made RFI with specific 
characteristics. The quiet period T prev rejects pulse trains. 
The secondary threshold T2, which is generally set lower 
than the primary threshold 7\, rejects signals with after- 
pulsing, even if the after-pulsing is weaker than the initial 
pulse. The condition on the timing and number of T2 
crossings allows for some ringing of a bandwidth-limited 
pulse but rejects long and/or irregular pulse trains or dig- 
ital spikes. The condition on Q = P ma x/NC can select on 
pulse shape, but in practice this is not needed when the 
other trigger parameters are set correctly. 

Typical values for the thresholds and noise rejection 
parameters are as follows: T\ = 200 ADC counts and T2 = 
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Figure 5: Sample LI triggers, showing the residual from a 20 ms 
period (50 Hz) vs. time. The drifting phase can be tracked with a 
digital PLL (black line). 

150 ADC counts, on an RMS noise level of 40 ADC counts; 
NC m in = 1 and NC max = 8 during a period of T per = 
6.25 /is; maximum crossing time TC max of 130 ns; and a 
quiet pre-trigger period of T prev = 1.25 /is. 

The algorithm can trigger on either signal polarization, 
or the logical AND or OR of the two. The narrowband fil- 
tering described in Sect. 13. H is performed before the trigger 
algorithm, but one can also trigger on the unfiltered signal; 
the unfiltered signal is stored in either case. 

3.3. Periodicity Filtering 

Man-made pulsed RFI, unlike cosmic-ray air showers, 
is often periodic in nature. In particular, by examining the 
time between successive level 1 triggers, clear signatures 
are visible at 100 Hz and sub-harmonics thereof (twice the 
power grid frequency of Argentina). 

One difficulty in vetoing such events is the drift in 
phase of the 100 Hz pulses. This can be tracked in the 
RDS software with the use of a digital phase-locked loop, 
which can track the changing phase and allow a veto of 
the periodic events (see Fig. [5j). 

3.4- Directional Filtering 

As shown in Fig. [2j most of the pulsed background 
events come from man-made sources on the horizon. When 
a level 3 trigger is formed with three or more stations in co- 
incidence, the DAQ has enough information to reconstruct 
the direction and veto these events, either with a simple 
selection on the reconstructed zenith angle, or by using 
additional azimuthal and timing information. The only 
limitation is the computation time required for a direc- 
tional reconstruction, which must keep up with the event 
rate. A plane- wave fit to the trigger times of the stations 
is sufficiently fast and accurate for such purposes. 

A faster method that can veto directional hot spots, 
but does not require a full directional reconstruction, uses 
the distribution of trigger time differences between two 
stations. Directional hot spots can be identified by peaks 
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Figure 6: Trigger-time difference between stations 1 and 4 (solid: 
sample data, showing a pulsed RFI point source; dashed: expected 
cosmic-ray distribution, normalized to the same number of events). 



in the time difference histograms for each station pair, and 
these distributions are quite different from that expected 
from air showers (see Fig. [6j). 

To remove the hot spots, the DAQ forms dynamic his- 
tograms F(i,j) of trigger time differences Ti — Tj for each 
station pair (i, j); the histograms continuously update us- 
ing a fixed number of past entries. For a candidate level 
3 trigger event, the number of "good" time differences 
(Ti , Tj ) is calculated as 
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if F(Ti — Tj) < F c 



otherwise 



(3) 



where F c is an adjustable threshold that depends on the 
bin size and number of entries in the histograms. The 
selection criterion to keep the event is N g00( ^ > N pa i rs /2 . 
Future revisions may include the expected time difference 
distribution of air-shower events as well. 

3.5. Offline Filtering 

To date, confirmed cosmic-ray events in AERA have 
been found by searching offline for time coincidences with 
the water- Cherenkov surface detector of the Pierre Auger 
Observatory. However, we are developing methods to iden- 
tify cosmic-ray candidates without information from other 
detectors. 

One promising technique is the use of the signal po- 
larization. Since the primary component of the air-shower 
radio emission is geomagnetic in nature [7], the direction 
of the predicted electric field e g is parallel to v x 6, where 
v is the direction of the shower axis, and b is the direction 
of the local magnetic field of the Earth. By deconvolving 
the response of the antenna and the station electronics, 
the measured electric field of a radio pulse E(t) can be 
reconstructed [8]. The direction of polarization can be 
compared to that expected from an air shower with purely 
geomagnetic emission, i.e., events are selected if 
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where Sfr c is a threshold space angle. This method does 
not yet account for the contribution of secondary emission 
mechanisms with different polarization signatures, such as 
that from the charge excess in the shower 0, [Toj]; however, 
this can in principle be added using additional information 
such as the core position of the shower. 

4. Future Plans 

AERA will expand in 2012 and 2013 to instrument 
an area of 20 km 2 with 160 radio detector stations. This 
expansion will entail a transition from communication over 
optical fiber to a wireless network. Several designs are 
under development for the wireless system, including: 

1. a fully custom time division multiple access (TDM A) 
system in the 2.4 GHz band that supports up to 180 
subscribers per channel, at a bandwidth of 5.5 Mbps; 

2. a commercial 802.1 In + TDM A system in the 5 GHz 
band that supports 80-100 subscribers per channel, 
at a peak bandwidth of 150 Mbps (80 Mbps typical 
TCP/IP throughput); and 

3. a distributed "gossip" protocol in which stations com- 
municate with each other and form multi-station co- 
incidences without communication with a central DAQ. 

The first two systems are currently being field tested at the 
AERA site, while the third system is in the development 
and simulation stage. 

In addition to the developments in data acquisition, 
triggering, filtering, and communications systems described 
here, the future stages of AERA will also use enhanced 
antennas, low- noise amplifiers, and digitization hardware. 
The increase in efficiency as well as detection area will 
result in a significant sample of hybrid and super-hybrid 
events that will be used to quantify the potential of the 
radio technique for measurement of cosmic-ray air shower 
energy, direction, and primary composition. 
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